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The predictions of a nonequilibrium schematic mode-coupling theory developed to describe the 
nonlinear rheology of soft glassy materials have been numerically tested in a sheared binary Lennard- 
Jones mixture. The theory gives an excellent description of the stress/temperature 'jamming' phase 
diagram of the system both at the microscopic and the macroscopic levels. In the present paper, 
we focus more particularly on the issue of an effective temperature T^s for the slow modes of the 
fluid, as defined from a generalized fluctuation-dissipation theorem. As predicted theoretically, many 
different observables are found to lead to the same value of T^s, suggesting several experimental 
procedures to measure Tcs- New, simple experimental protocols to access T^b from a generalized 
equipartition theorem are also proposed, and one such experiment is numerically performed. These 
results give strong support to the thermodynamic interpretation of Tes and make it experimentally 
accessible in a very direct way. 
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Glassy materials are usually defined by the fact that 
their relaxation time is larger than the experimental time 
scale. In simple molecular systems, the associated glass 
transition temperature corresponds to very high viscosi- 
ties, making it difhcult to investigate experimentally their 
rheological properties. In complex fluids (e.g. colloids, 
emulsions) it is however possible to reach a glassy situa- 
tion, in the sense of large relaxation times, with systems 
having viscosities or shear moduli that allow for rheologi- 
cal investigations |l| . Such materials have been described 
as 'soft glassy materials' |Q. 

In its glassy state, a material is by definition out of 
equilibrium. Physical properties are then a function of 
the time t^, spent in the glassy phase, a behavior called 
aging Interestingly, recent experiments on various 
complex fiuids have demonstrated striking similarities 
with other, more standard, glassy systems Upon 
imposing a steady, homogeneous shear flow, a different 
kind of nonequilibrium situation is obtained. The flow, 
characterized by the shear rate 7, creates a nonequilib- 
rium steady state, in which time translation invariance 
is recovered jo). This situation can therefore be used to 
probe the glassy state, with the convenient feature of hav- 
ing the shear rate 7 rather than the waiting time tw as 
a control parameter. Moreover, this way of probing the 
nonequilibrium properties of glassy systems is probably 
more relevant experimentally than the aging approach, 
at least in the case of soft glassy materials. 

Recently, a general scenario was proposed for glassy 
systems subject to an external forcing based on the 
study of mean-fleld models. The rationale of this ap- 
proach is that the equilibrium dynamics of these mod- 
els is equivalent to the 'schematic' mode-coupling ap- 
proach of slowing down in supercooled liquids pT| , p^ . 



The study of their nonequilibrium dynamics can thus be 
seen as a nonequilibrium schematic mode-coupling ap- 
proach |H,|l3|,0. To our knowledge, the mode-coupling 
theory of supercooled fluid has not yet been extended 
to fluid under shear beyond the linear response regime. 
However, in analogy to what was done for aging or su- 
percooled systems, it is sensible to bypass this aspect 
and to carry out a direct comparison between mean-field 
predictions and experimental or numerical results. 

The aim of this work is thus to check on a realistic 
model of a fiuid the predictions that emerged from the 
theoretical approach of Ref. Several earlier stud- 

ies have been devoted to simulating glassy sheared flu- 
ids. Ref. focused on the dynamics at the molecular 
level. Our study is devoted to more global aspects, with 
the aim of providing experimentally testable predictions. 
Ref. [|6| also proposed to use the shear rate as a control 
parameter for jamming systems. Some dynamic aspects 
of a model similar to ours, with the difference of being 
athermal (zero temperature) have been investigated, as 
discussed in a companion paper . 

We have investigated the stress/temperature jamming 
phase diagram (Fig. |l|) of the 'standard model' for super- 
cooled liquids, namely a 3D 80:20 binary Lennard- Jones 
mixture, in a simple shear ffow defined by v = jyex- The 
system consists of = 2916 particles in a cubic simu- 
lation box. It has been characterized in much details at 
the reduced density p — 1.2, where we carry out our sim- 
ulations, both in supercooled [Q and aging ||l^ regimes. 
Our simulations follow the protocol detailed in Ref. : 
the system is first made stationary on a timescale of a 
few 7^^. Then we perform our measurements in a range 
T e [0.15,0.6] and 7 G [10-^,10-1]. Standard Lennard- 
Jones units |jl^,^ are used. 
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Following the theoretical analysis of Ref. [T^] , our ap- 
proach consists in three main steps. We first investigate 
macroscopic rheological aspects (flow curves). We then 
focus on microscopic properties by an analysis of the den- 
sity fluctuations: structure factor, density-density corre- 
lation functions. Both aspects will be discussed in much 
more detail in a longer publication Last, we anal- 

yse the main new feature for the field of rheology which 
emerged in Ref. namely the existence, behavior and 
properties of an effective temperature, defined through 
a generalized fluctuation-dissipation theorem p^ . This 
point, together with its experimental consequences are 
the main object of this paper. 




FIG. 1. The {o,T) plane of the jamming phase diagram. 
The dashed curves are the viscosity contour plots with 77 = 20, 
30, 50, 100, 200, 500, 1000 and 2000 (from top to bottom). 
The full line is the yield stress cro(T), from Eq. (|l]). Arrow 
marks the mode-coupling temperature Tc ~ 0.435. 



The macroscopic behavior of the system under shear 
is summarized in a stress/temperature phase diagram in 
Fig. |l|. Our results are qualitatively similar to the one 
reported in Ref. jl^ on a different glass-former. A New- 
tonian behavior is observed in the high-T, low-cr part of 
the phase diagram. This corresponds roughly to situ- 
ations where 7"^ is larger than the relaxation time of 
the fluid, which is thus only weakly affected by the flow. 
Outside this region, we find that the viscosity of the fluid 
decreases when the shear rate increases, a shear-thinning 
behavior reported in various 'soft glassy materials' [Q. 
We used the phenomenological relation [n| 



(To + ai 



(1) 



to extract the behavior of the yield stress ctq (T) , reported 
in Fig. ^ This 'jamming transition' line was recently 
experimentally investigated ||2^. Although our results 
cover 3 decades in shear rate, we cannot report a defini- 
tive functional form for the shear-thinning of the vis- 
cosity [p3| which can also be satisfactorily described by 



CTo = 0, which amounts to describe the system as a power- 
law fluid, rj ^ 7""^. In the supercooled regime we find 
n ~ 1/3, while at lower temperatures, the exponent n is 
temperature dependent with n ^ when T — > 0, indicat- 
ing that a yield stress could exist in this limit only 
The latter power-law behavior is precisely the one pre- 
dicted theoretically in Ref. 



|10| 



Equilibrium mode-coupling theory gives rise to several 
quantitative predictions regarding the scaling properties 
of the intermediate scattering function, 



Ck{t) 



1 ^ 

- ^ (exp {ik ■ [rj{t + to) - rjito)])) , (2) 



when the glassy phase is approached by lowering T with 
(7 = 0. It was shown in Ref. ||l^] that similar scal- 
ing properties are expected when the glassy phase is 
approached by lowering the shear stress cr at constant 
temperature (vertical line in Fig. |l|). We have tested in 
a detailed way these predictions in our simulation 
Here, we only report the validity of the theoretically pre- 
dicted 'time-shear superposition property' jl^ in Fig. ||, 
which shows that the slow decay of Ck{t) has the scal- 
ing form Ck{t) ~ F(t/t^f,i): where the relaxation time is 
defined, as usual, by Cfc(trci) = e^^- The scaling func- 
tion F{x) is well described by a stretched exponential, 



F{x) ~ 
from f3 ' 



exp(— a;f^), with an exponent /3 which increases 
- 0.77 for T > to the value /3 = 1 as T 0. 
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FIG. 2. Correlation functions for T = 0.3 < Tc and differ- 
ent values of the shear rate, 7 = IQ-*, 3 ■ 10-'', 10-^, 3 • IQ-^, 
10-^ 3•10-^ 10-^ (from left to right at C = 0.9) can be col- 
lapsed if the time is rescaled by trci(7). The dashed line is a fit 
to a stretched exponential form, with an exponent /3 = 0.95. 



We now focus on the issue of an effective temperature 
for the slow modes of the sheared fluid. This quantity, 
naturally included in nonequilibrium mode-coupling the- 
ories pO| , is defined through a nonequilibrium generaliza- 
tion of the fluctuation-dissipation theorem (FDT) O]. 
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Consider two physical observables, 0{t) and 0'{t), their 
connected cross-correlation function Coo'{i) = {0{t + 
to)0'{to)) — {0{to)){0'{to)), and the response function 
Roo'it) = ^sh^^i\to)'^ ' where ho' is the field thermo- 
dynamically conjugated to 0'{t). At equilibrium, both 
quantities satisfy the FDT, TRoo'{t) 



dCr 



'it) 



dt 



The 

susceptibility xoo' (i) = /J dt' Roo' {t') can be measured 
by applying a small (to ensure linear response), constant 
field ho' between times and t and FDT implies a simple 
linear relation Txoo'{t) — (CoO'(O) ~ C'oo'(^))- the 
sheared fluid, an effective temperature is defined by 



Roo'{t) — 



(3) 



where the T^'^ (x) are a priori arbitrary functions of 
their argument, which may in general depend on the ob- 
servables O and O' under study. For any pair of observ- 
ables, (x) can be measured in the sheared system, 
by following the same linear response procedure as above, 

so that XOO' (t) = /( 



'(0) 



Coo' it) TOO' (^) 



The existence of an 



effective temperature is thus demonstrated if a straight 
line is obtained in a susceptibility-correlation plot pa- 
rameterized by the time. Obviously, the introduction of 
an effective temperature is of 'thermodynamic' interest 
only if this quantity is actually independent of the ob- 
servables under consideration, = Tcff- This crucial 
feature is true at the mean- field level | p^ , p^ , and we shall 
prove that it is nicely satisfied in our model. 

We have already shown the existence of such an effec- 
tive temperature for the slow modes of the system for a 
single correlation function, Eq. (|^), at a given wavevec- 
tor We also found that the essential phcnomeno- 

logical idea that a system sheared more vigorously has a 
higher Tcff is indeed captured by this definition |2C| ]. 

Here, we go much further in our investigations and 
prove that several different observables lead to the same 
value of Tefi . Since the numerical measurement of Teff is 
very demanding, our strategy has been to compute for 
T = 0.3 and 7 = IQ-^ the value of Teff from the ob- 
servables of Ref. with a great accuracy. We found 
Teff ~ 0.65. We then computed Tes using different ob- 
servables, and checked that the value so obtained was 
compatible with Tcff — 0.65. Our results are summa- 
rized in Fig. ^ which shows 10 among the 14 differ- 
ent susceptibility-correlation measurements we have per- 
formed. All our data are well compatible with the single 
value of Teff = 0.65 for the slow modes of the fluid. 

We first investigated the density fluctuations as an ob- 
servable, taking 0{t) = j; ^jLi ^ j {^^ ■ ''^(t)); 
0'{t) = 2^^-^ £j cos (fc • rj(t)). For some wavevectors, 
this was done separately for both types [A and B) of 
Lennard-Jones particles (noted 'A + S' in Fig. ||). Taking 
Ej = 1 selects the coherent part of the intermediate scat- 
tering function (noted 'coh' in Fig. while £j = ±1 se- 



lects the incoherent one. Such response-correlation plots 
could be obtained experimentally in soft condensed mat- 
ter systems. While the correlation functions are reason- 
ably easily obtainable through light scattering experi- 
ments, the same is not true of response functions. To 
obtain the latter, one has to manipulate the particles 
through some externally applied potential, modulated at 
the same wavevector as used in the light scattering exper- 
iment. One suggestion would be to use some non- index 
matched tracer particles in an index matched colloidal 
suspension. The tracer particles would then be sensitive 
to the intensity of the local electric field, as in optical 
tweezers. An interference pattern would actually realize 
the modulated external potential considered here. Read- 
ing of the response could then be obtained from the scat- 
tering at the wavelength corresponding to this pattern. 
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FIG. 3. Ten representative susceptibility-correlation plots 
are shown to be consistent with the same Teff for the slow 
modes (small value of the correlation). Full lines are the equi- 
librium FDT of slope —l/T, dashed lines have slope —l/Tefs 
with Tcff — 0.65. Numbers refer to wavevectors (see also the 
text). 

We also used as a correlation the mean square dis- 
placement of a tagged particle. The associated response 
function is the displacement induced by applying a small, 
constant external force to this tagged particle . Both 
quantities are linked by a FDT, the Einstein relation. 
We computed both quantities separately for particles of 
type A and B (noted 'Self A^ and 'Self B' respectively). 
Its interest, especially in view of experimental realiza- 
tions, lies in the fact that the full time dependence of 
the correlation and response functions is not needed to 
extract Toff. Indeed, at large times, both quantities be- 
come proportional to the time which defines the diffusion 
constant and the mobility. One may therefore define Toff 
simply as the ratio of diffusion to mobility. Again, exper- 
iments could be considered if tracer particles sensitive to 
an external force field (e.g. magnetic particles) could be 
introduced into the system, and their mobility measured 
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together with their diffusion constant. 

A completely different observable, relevant to flow situ- 
ations, is the stress a. We have studied the case in which 
the observables O and O' are equal to the diagonal stress 
in the direction transverse to the flow, Uzz{t)- To add a 
field conjugated to Uzzi a- compression 8Lz of the simula- 
tion box is realized by rescaling all particle coordinates at 
time i = 0, in analogy to stress relaxation experiments. 
The corresponding curve is labelled 'Stress' in Fig. |. 
Experimentally, the off-diagonal component of the stress 
would be used as the observable. Preliminary results in 
this direction, using an extremely sensitive rheometer, 
have been obtained in aging colloidal systems [^ . 
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FIG. 4. Mass dependence of the mean kinetic energy in the 
z direction for T = 0.3 and 7 = 10"^. Horizontal lines are 
T = 0.3 and Tcff = 0.65. Error bars are evaluated from tracer 
to tracer fluctuations. 



fluid, e.g. polymeric, or by investigating rotational de- 
grees of freedom of non-spherical tracer particles. 

This last result opens the way for new, simple deter- 
minations of Toff in out of equilibrium glassy materials. 
Indeed, no 'complex' dynamic functions such as correla- 
tions or susceptibilities are needed here. In this context, 
it would be interesting to reproduce Perrin's experiment 
on barometric equilibrium of colloidal suspensions [ p6[ . 
We expect indeed that the barometric equilibrium of 
heavy particles inside a horizontally sheared fluid should 
be ruled by T^s instead of the room temperature. 

All these results give strong support to the theoreti- 
cal scenario elaborated from mean-field theories to de- 
scribe the rheology of soft glassy materials jl^ . Note in 
particular that Teff, which we have shown to be a phys- 
ically relevant quantity, is a natural outcome of the the- 
ory. While numerical simulations provide a test of the 
theory on short time scales, systematic experiments that 
could test quantitatively existing theories of nonequilib- 
rium glassy dynamics are still needed. In that sense, 
the situation is more or less similar to the mid-eighties 
when schematic equilibrium mode-coupling were already 
derived, but with little experimental confirmation of its 
main features. This is why we tried, as much as possible, 
to suggest experimental counterparts to our numerical 
measurements. We hope that our findings and sugges- 
tions will motivate further experimental work in the field. 

We thank L. Bocquet, J. Kurchan and W. Kob for 
discussions. This work was supported by the PSMN at 
ENS Lyon and the CDCSP at Universite de Lyon. 



Definition (||) of Tcff implies that, if the fluid is used as 
a thermal bath to equilibrate a subsystem (a 'thermome- 
ter') of typical time scale ^ <i.ob the thermometer does 
not measure the microscopic temperature, but rather the 
effective temperature associated with its characteristic 
time scale We propose here to use tracers of mass 
mtr as a thermometer, since tuning mtr allows to con- 
trol their vibration time scale tg. We have considered 10 
massive tracers with mtr G [IjlO^], but being otherwise 
identical to A particles. Since ts ~ ^mtr, heaviest parti- 
cles have a frequency typically 10^ times smaller than the 
light ones, implying ta ~ t-cei- Reading of the temperature 
is done by measuring the average mean square velocity 
of the tracers in the direction z. Results are shown in 
Fig. ^, which shows that light particles measure the bath 
temperature while heaviest particles measure Teff. This 
implies that a generalized equipartition theorem holds, 

(^\mtrvl^ = ^Teff. (4) 

This result could be tested against experiments involving 
for instance colloidal particles, the tracers, in a complex 
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